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ABSTRACT 



O . 

We investigate the physical and chemical conditions necessary for low-mass 
O ' star formation in extragalactic environments by calculating various character- 

istic timescales associated with star formation for a range of initial conditions. 
The balance of these timescales indicates whether low-mass star formation is en- 
Oh 1 hanced or inhibited under certain physical conditions. In this study, we consider 

timescales for free-fall, cooling, freeze-out, desorption, chemistry and ambipolar 
■ diffusion and their variations with changes in the gas density, metallicity, cosmic 

ray ionisation rate and FUV radiation field strength. We find that extragalac- 
tic systems with high FUV radiation field strengths and high cosmic ray fluxes 
considered at a range of metallicities, are likely to have enhanced low-mass star 
formation unless the magnetic pressure is sufficient to halt collapse. Our results 
\Q • indicate that this is only likely to be the case for high-redshift galaxies approach- 

ing solar metallicities. Unless this is true for all high-redshift sources, this study 
finds little evidence for a high-mass biased IMF at high redshifts. 



Subject headings: astrochemistry — ISM — stars: formation 



1. Introduction 



Recent multi- wavelength studies of galaxies have permitted es timates to be made of the 



evolution of the space density of the massive star formation rate (IHopkins fc Beacoml 12006 



Lilly et all 119961 ; iMadau et al.lll996l ). While these rates are at present poorly constrained, 
there is currently no information on the formation at high redshift of low mass stars, a process 
that in the Milky Way is known to accompany the formation of massive stars. The low mass 
stars may be optically insiginificant, but they modify the physics of the host galaxies rather 
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significantly. In particular, the shape of the stellar initial mass function (IMF) and its slope 
as specified by the ratio of low-mass to high-mass stars, are crucial for determining the 
star formation history and stellar mass dens i ty of a galaxy. The r e has been much interest 



in va rying the stellar IMF (jEhriegreenl 120081 ; IWilkins et al.l 120081 ; iDavel 120081 ; Ivan Dokkum 



20081 ) and recent studies have suggested that the IMF must be top-heavy at high redshifts 
to explain the counts of faint subm illimeter galaxies and the luminosity function of Lyman 
break galaxies (IBaugh et al.ll2005l ). However, this has by no means been proven true for all 
high redshift environments. 

Studies of low mass star formation in the Milky Way have shown that it arises in the 
interplay of a number of competing processes including gravitational collapse and thermal 
and magnetic support. Thermal support at the low temperatures of interstellar clouds is 
determined by the cooling mechanisms, while magnetic support is controlled by the ionization 
fraction; both are largely determined by chemistry which itself is moderated by cosmic ray 
and photon fluxes, and metallicity. A consideration of the timescales appropriate for the 
Milky Way indicates - as it should - that low mass star formation takes place under Milky Way 
conditions. However, the physical properties in high redshift galaxies may depart significantly 
from the Milky Way values. In this paper, therefore, we compute in a rather crude way how 
these relevant timescales may vary when these physical conditions are modified from Milky 
Way values. We use the results to determine the regions of parameter space in which low 
mass star formation is unlikely to occur. By doing so, we are able to investigate claims of a 
high-mass biased IMF at high redshifts. 

In § [2] we list the relevant physical processes and give estimates of the timescales. The 
cooling functions and ionization fractions are determined from a complex model of interstellar 
chemistry described in § [31 Our results are described in § HJ and some general conclusions 
regarding low mass star formation at high redshifts are made in § 



2. Timescales for Star Formation 

Star formation is a complex process. While gravity is the main driver, there are many 
other mechanisms that compete to aid and oppose the collapse of molecular clouds. Cooling 
by molecular radiation leads to thermally induced collapse of a molecular cloud and if the 
timescale associated with this process is longer than the free-fall time, the collapse will be 
halted. The rate of cooling in turn is determined by the abundances of the various different 
molecular coolants such as CO and so we can see the role that chemistry plays. The chemistry 
itself is driven by cosmic ray ionisation of molecular hydrogen. In some circumstances, the 
surfaces of dust grains can become coated with ice due to the freeze out of different species 
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onto the grains. This process is effective in removing coolant molecules from the gas phase. 
Energy released from H 2 formation can liberate the frozen species into their gaseous phase. 
Interstellar magnetic fields too can affect the collapse of molecular clouds. The ions are 
coupled to field lines and neutral atoms to the ions via collisions. The fractional ionisation 
therefore determines the magnetic pressure and the level of ion-neutral friction which in turn 
can halt collapse. 

In the local universe, the timescales associated with these var ious competing processes 



are found to be roughly of the same order of magnitude - ~ 10 5 yrs (jWilliams!ll998l ) for a core 
density of ~ 10 4 cm -3 . This balance of timescales is clearly conducive to star formation. By 
exploring the change in these timescales for a range of initial conditions, we can effectively 
explore different regions of redshift space and hope to learn something about when and how 
the first low mass stars started to form. 



2.1. Collapse 

Low mass star formation essentially occurs by the collapse of a molecular cloud under 
its own weight accompanied by the dissipation of gravitational potential energy. In any 
situation where the star formation is not externally triggered by a violent event, for example 
the passage of a shock, the actual collapse time will be somewhat longer than the free-fall 
time due to turbulence within the cloud. In this work, we use the free-fall time as a useful 
under-estimate of the typical collapse time for a molecular cloud undergoing quiescent star 
formation. 

We can relate the free-fall time of a collapsing core to the number density of hydrogen 
nuclei as follows: 



37T 1 

'-ff = \/ 777777 = °- 75 x W 8 /{n H T 2 yr (1) 



16Gp 

where = (n(H) + 2n(H 2 )) is the local number density of hydrogen nuclei in cm 



-3 



flDvson & Williamsl Il997f ). 



2.2. Cooling 



In order for a cloud of gas to continue to collapse and form a star, the cloud must continue 
to lose energy, otherwise, the high temperatures and densities within the cloud would halt 
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collapse before the formation of a star could occur. Energy loss occurs via different atomic 
and molecular transitions, while at high densities, collisions between hot gas and cooler dust 
grains can also be efficient in cooling the gas. The rate at which energy is lost depends on 
the cooling function, A tot / erg cm ~ 3 s^ 1 , which in turn depends on the abundance of various 
different molecular coolants. The total cooling time is then given by: 



tcool = ~rr yr (2) 

where k is the Boltzmann constant, T is the gas temperature and £ is the metallicity in 
units of solar metallicity. The quantity, |n^/cT represents the total energy, both thermal 
and internal, of the system and £A to ( is the total interstellar cooling rate. 

Note that our chemical models described in § [3j do not include the effects of remov- 
ing molecular coolants from the gas phase by freeze-out. If freeze-out occurs, the cooling 
timescale will be underestimated in regions of parameter space where there is significant 
ice-formation. 



2.3. Freeze-out 



A d etailed calculation of th e rate per unit volume at which a species freezes out can be 
found in iRawlings et all (1l992f ) . This quantity depends on the relative abundance of dust 
grains compared to hydrogen nuclei, the physical properties of the grain as well as the species 
whose freeze-out is being considered, and electrostatic effects to take into account the fact 
that ionic species will freeze-out at a different rate to neutral species. Due to the presence of 
more than one grain population in most systems , the g rain properties need to be averaged 
over the grain size distribution. IRawlings et al.l (119921 ) relate this size distribution to the 
depletion coefficient, D which depends on the abundance of very small grains and on the 
metallicity. 



Taking these assumptions into consideration, t 



re freeze-out timescale can be derived 



from the rate of freeze-out given by IRawlings et al.l (Il992l ) : 



l fo 



0.9 x 10 6 (m x /28)2 
(T/10)*(n H /10 4 )D 



yr 



(3) 



where mx is the molecular mass of the species whose freeze-out is being considered and 
T is the gas temperature. 
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2.4. Desorption 

Desorption is the process by which molecules in ices on the surface of dust grains, are 
returned to the ISM. The mechanism through which heat is generated for this desorption, 
is still debated ([Roberts et al.l 120071 ). In this work, we consider only one mechanism for 
desorption, namely exothermic reactions occurring on the grain surface. More specifically, 
the liberation of energy from H2 formation on dust grains results in a proportion of the 
frozen molecules on the dust grain being returned to the gas phase. 

We calculate the desorption rate of CO which is the most abundant molecular species 
besides hydrogen. This rate is: 



The rate depends on the metallicity which scales linearly with the dust to gas ratio, 
the abundance of both molecular and atomic hydrogen, the temperature, T and the number 
of CO molecules desorbed for every H2 formed, 7. The coefficient 3 x 10~ 17 represents the 
canonical H2 formation rate in the Milky Way. 

From this expression, we can calculate a typical timescale associated with desorption: 

1 x 10 9 



l des 



[n H n(H)(T/W0) l n ^] 



yr (5) 



where n# is the local number density of hydrogen nuclei and n(H) is the column density 
of hydrogen atoms only. The number of CO molecules desorbed per H 2 formation, 7 can be 
constrained by considering a pre-stellar core in the Milk y Way like L1689B w here significant 



depletion of CO molecules has been directly observed (IRedman et al.ll2002l ). Assuming for 
such an object with n# ~ 10 4 — 10 5 cm -3 , that the timescale for freeze-out is less than half 
the timescale for desorption gives an upper limit on 7 of ~0.25. We use 7 = 0.2 throughout 
this work and assume that this number does not change in high redshift systems. This is 
plausible as we expect the microscopic properties of the dust in high redshift galaxies to 
be fairly similar to those in the Milky Way. Both the freeze-out and desorption timescales 
depend on these microscopic properties. 
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2.5. Ambipolar Diffusion 

Magnetic fields may play an important role in supporting dense molecular clouds and 
many of the fragments within them. The magnetic field lines are coupled directly to the 
ions and indirectly to the neutral matter through frequent ion-neutral collisions. At low 
ionisations, the frequency of the ion-neutral collisions becomes less regular and the field 
lines decouple from the neutral matter thereby allowing the matter to drift in response 
to the gravitational potential. This process is known as ambipolar diffusion and plays an 
important role in determining the rate of quiescent star formation. 

The characteristic timescale associated with this process of decoupling the field from the 
neutral matter will depend o n the relative drift velocities of the ions and neutral particles. 



Hartquist fc Williams! (119891 ) define the ambipolar diffusion timescale as follows: 



^amb 



4 x 10 5 (xi/10- 8 )?/r 



(6) 



This is a function only of the fr actional ionisation, X j within the cloud and is consistent 
with the more detailed treatment of iMouschoviasI fll979l ). However, in the earlier paper, the 
ambipolar diffusion timescale is shown to be a function of position within the cloud. There- 
fore variations in the degree of ionisation as a function of cloud depth become important. 
For the purposes of our rather cr ude calculati o ns, w e neglect this positional dependence. 
Eq. [6] therefore corresponds to the IMouschoviasI (119791 ) timescale at an intermediate position 
towards the interior of the core. It should be noted that this timescale is only relevant for 
magnetically sub-critical cores. 



2.6. Ion-Molecule Chemistry 

The chemical timescale is the time taken to create molecular coolants from the initially 
atomic gas assuming that nearly all of the hydrogen is molecular. In other words, it is the 
time taken to ionise sufficient hydrogen to allow complete chemical conversion of C and 
O to molecular form. Cosmic rays are responsible for ionising the hydrogen. The relative 
abundances of carbon and oxygen will scale with the metallicity, £. The chemical timescale 
is 



3£(ng* + ng*) 5 x 10 6 £ 
n(H 2 )C ~ C/l x 10- 



lchem ~ ~(Ti\f~ ^/1 vm-iT^ ( 7 ) 



The chemistry is important for providing the coolant molecules and therefore must not 
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be too long if star formation is to occur. The analytical expression for the chemical timescale 
represented by Eq. [7| depends only on the metallicity and cosmic ray ionisation rate and is 
rather crude. The chemical timescale is less relevant in cores than in the parent cloud which 
may be collapsing from a diffuse to a dark state as there isn't much change in the chemistry 
once cores begin to form. 



3. Modelling 

In order to model the various different timescales discussed in § [2] we use a time and 
depth-dependent chemical code for molecular clouds and photon-dominated regions to solve 
simultaneously for the chemistry, thermal balance and radiative transfer within a cloud. 
The code models the gas phase chemistry inside the cloud which is treated as a ID semi- 
infinite slab o f gas illuminated by a uni-directional flux of FUV photons incident on its 



outer surface (jBell et al.l 120061 ). Differences in model geometry are unlikely to affect the 



chemi stry significantly enough to influence our rather qualitative results (See iRollig et al. 



(120071 ) for a PDR benchmarking code com parison) . Th e FUV radiation field is taken to 



be the standard interstellar radiation field (jDrainelll978l ) multiplied by a factor of % where 
X = 1.7 corresponds to a flux of 1.6 x 10" 3 erg s -1 cm -2 , appropriate for the unshielded 
interstellar medium in the Milky Way. The radiation field strength declines exponentially 
with increasing depth within the cloud and therefore depends on the visual extinction A v . No 
grain chemistry is included and hence the freeze-out and desorption of species is ignored in 
our models. Various different heating and cooling mechanisms are invoked within the code. 
Photoejection of electrons from dust grains and polycyclic aromatic hydrocarbons (PAHs) 
and FUV pumping and photo-dissociation of H 2 are the dominant heating mechanisms while 
the gas is mainly cooled through emission from collisionally excited atoms and molecules and 
by interactions with the cooler dust grains. The gas temperature is determined at each depth 
point by using an iterative procedure to balance the total heating and cooling rates. Although 
the chemical model is also time dependent, we consider all outputs at steady state, assumed 
to correspond to a time of lGyr. 

The code takes as its inputs the number density of hydrogen nuclei, n#, the metallicity, 
£, the cosmic ray ionisation rate, (, the strength of the incident FUV radiation field, \ and 
the turbulent velocity, v within the cloud. As we are considering quiescent low-mass star 
formation only, we fix the turbulent velocity at a value of v = 1.5 km s" 1 . We will consider 
variations in the other parameters and implications for the star formation timescales in § HJ 



The chemistry consists of a chemical network containing 83 species. Thes e species inter- 
act in over 1300 reactions, whose rates are taken from the UMIST99 database (]Le Teuff et al. 
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2000l ). Initially, hydrogen is assumed to be in atomic form. 



Given these inputs, the model calculates chemical abundances for the different species, 
emission line strengths and gas temperatures as a function of the visual extinction, A v for 
every time step. 



4. Results 

In this section we describe how the various timescales associated with low mass star 
formation described in detail in § [2] vary as we change the free parameters. We consider 
variations in the number density of hydrogen nuclei, n#, the metallicity, £, the cosmic ray 
ionisation rate, £ and the strength of the incident FUV radiation field, \. When varying 
any of the four parameters, the others are fixed at their local values namely nn = 10 5 cm -3 , 
£ = £©; C = s" 1 and x — 1-7- All outputs are considered at steady state corresponding 

to a time of lGyr and at two depths corresponding to A v ~ 3 and A v ~ 10. Outputs of the 
code include the cooling function, A tot , the fractional ion abundance, Xi, the gas temperature 
T and the column density of hydrogen atoms n(H) and these values can be substituted into 
equations[2],[3l[5]and[6]to calculate the cooling, freeze-out, desorption and ambipolar diffusion 
timescales. 

For each set of initial conditions, we make some references to extragalactic objects that 
may exhibit these properties and hence are able to infer if low mass star formation is likely 
to be encouraged or inhibited in such objects. We emphasise that our calculations are crude 
in nature and instructive only in terms of observing general trends. 



4.1. Varying the Gas Density 

We vary the hydrogen nuclei number density between 100 cm -3 and 10 6 cm -3 in this 
study and look at the effect of this variation on the low mass star formation timescales. The 
results are summarised in Figure [lj 

As expected, the free-fall time decreases with increasing number density converging to 
~ 10 6 yrs at 10 4 cm" 3 . At this density, all the different timescales converge to roughly the 
same order of magnitude. Thus both freeze-out and desorption proceed with desorption 
a slightly slower process so that molecular ices do form. Desorption maintains a popula- 
tion of molecular coolants in the gas phase. The free-fall time is comparable to the cooling 
time, thermally induced collapse occurs and stars are able to form. The ambipolar diffusion 
timescale is comparable to the dynamical free-fall time and star-formation can occur in mag- 
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netically sub-critical couds. We note that these conditions are appropriate for a Milky Way 
type galaxy and therefore our approach indicates as expected that low mass star formation 
is reasonably efficient in Milky Way type systems. 

However, at very high densities, the free-fall time is shorter than the cooling time and 
the cores will heat up during the collapse which is then arrested. 

The freeze-out timescale and desorption timescale are comparable and both decrease 
with increasing density at the same rate suggesting that molecular ices are able to form 
at all densities between 100 cm -3 and 10 6 cm -3 . This also implies that the cooling time 
is systematically under-estimated as ice formation, which is not included in our chemical 
model, removes molecular coolants from the gas phase at all densities. 

Looking deeper into the cloud, i.e. going from A v ~ 3 to A v ~ 10 does not seem to affect 
significantly the freeze-out or desorption timescales. However, the cooling time is slightly 
larger at higher visual extinctions as the cooling function decreases with increasing A v . This 
is probably due to increased radiation trapping in the cooling lines. 

The ambipolar diffusion timescale is much longer than the other timescales at very 
low densities. This suggests as expected that star formation cannot occur in clouds with 
Uh < 10 3 cm -3 that are magnetically sub-critical as the field is very slow to decouple. 

4.2. Varying the Metallicity 

The metallicity, £, is another important physical property that will affect low mass star 
formation. This parameter is particularly important for modelling high redshift systems as 
the primordial gas from which these objects formed is likely to be less chemically enriched 
than the interstellar medium in the Milky Way in some cases, and will therefore have lower 
metallicities. As we are interested in low metallicity objects that formed early in the history 
of the Universe, we consider metallicities as low as 1O _4 £ . We also consider metallicities 
enhanced above solar values to 3£ . In our model, we assume that the elemental abundances 
of all metals scale linearly with metallicity. We also assume that the dust to gas mass ratio 
scales linearly with metallicity. The results from this study are summarised in Figure [2j 

As expected, the cooling time falls with increasing metallicity due to the increased 
abundance of molecular coolants. We can see that at metallicities below £ = O.l£ ,the 
cooling time is much longer than the free-fall time and collapse and star-formation are 
likely to be impeded by gas pressure. The freeze-out time also decreases with increasing 
metallicity due to the presence of more dust grains but the desorption time is independent 
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of metallicity. At low metallicities, freeze-out is slower than desorption so ice-formation is 
inhibited. Temperatures however remain high due to the long cooling timescales. 

The ambipolar diffusion timescale remains fairly insensitive to the metallicity. The 
fractional ionisation is much higher at low visual extinctions so the ambipolar diffusion 
timescale falls with increasing depth into the cloud. Note that the visual extinction scales 
linearly with metallicity so high visual extinctions in low metallicity clouds will correspond to 
larger physical distances within the cloud. Hence, the calculations for A v ~ 10 are probably 
unfeasible for low metallicity clouds. We have considered outputs at A v ~ 7.5 and A v ~ 9.5 
for £ = 1O _4 £ and £ = 1O _3 £ respectively. As the cloud is likely to be optically thick even 
at these lower depths, the chemistry is unlikely to be very different to the A v ~ 10 case. 

The chemical timescale increases linearly with metallicity as expected. Low metallicity 
systems therefore convert a significant fraction of available elements into molecular coolants. 



4.3. Varying the Cosmic Ray Ionisation Rate 



Although little is known about the origin of extragalactic cosmic rays, they are very 
important in ionising hydrogen and driving the chemistry during star formation. The level 
of ionisation within a dark cloud may determine whether or not the cloud is stable against 
gravitational collapse if magnetic fields are present. 



An estimate for the cosmic ray ionisation rate in the Galaxy is ~ 10~ 17 s" 1 flHartquist et al 



19781 ). In this study, we vary the cosmic ray ionisation rate between 10 18 s 1 and 10 13 s 1 



in order to mimic both quiescent star forming systems such as those in our own galaxy as 
well as starburst regions and AGN which are expected to have cosmi c ray fluxes of the orde r 
of 10 _15 s -1 to lO^V 1 jMeiierink et al.l bood : ISuchkov et all Il993l : Bradford et"aD l2Q03h . 
The results of varying the cosmic ray ionisation rate are summarised in Figure [3j 

As the cosmic ray ionisation rate is increased, the cooling time falls and at ( ~ 10~ 16 
s -1 , it drops below the free-fall time. Cores cool more and more quickly and are therefore 
able to collapse to form stars. The chemical timescale is small at high values of the cosmic 
ray ionisation rate explaining the low cooling times as molecular coolants are quick to form 
in such circumstances. 

The freeze-out time is greater than the desorption time and the propensity for ice- 
formation is reduced thereby allowing molecular coolants to remain in the gas phase to 
effectively cool the collapsing gas. Desorption is therefore unlikely to be a significant process. 
However, formally, the desorption time becomes shorter as expected with increasing cosmic 
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ray ionisation rate as more and more hydrogen molecules are ionised thereby increasing the 
abundance of atomic hydrogen as well as the temperature of the gas. If ices could be formed, 
they would therefore be quick to evaporate off the dust grains once again replenishing the 
gas phase with molecular coolants and encouraging star formation. 

The ambipolar diffusion timescale becomes very large as the fractional ionisation in- 
creases a nd the magnetic field is slo w to be release d from the collapsing core. Based on the 
results of ICrutcher et all (119871 ) and lAthreya et al.1 (119981 ). the magnetic field strength, B of 
molecular clouds varies between 1/iGauss and 100/iGauss for densities between 1cm -3 and 
10 5 cm" 3 . This fiel d strength, B scales with the density as B oc {nu) K where g < k < ^ 
( Mouschovias 1976b is consistent with observations . In this study, based on various ob- 
servations jCrutcher et al.l Il987l : iBourke et all hoOll : beck et al.l boOfit kirn et all Il99pl ) we 
consider k ~ 0.4. We can calculate the ratio of the magnetic pressure and thermal pressure, 
a using Eq. [HJ 



B 2 



a 



SnnnkT 



(8) 



For a » 1 the magnetic pressure dominates and for a << 1, the thermal pressure domi- 
nates. For high cosmic ray fluxes and at a constant density of n#=10 5 cm -3 , we compute 
a ~ 0.4. This suggests that the magnetic pressure may in some cases be sufficient to halt the 
collapse but in other cases, depending on the temperature of the core and the exact value of 
k, low-mass star formation can proceed unimpeded. 



4.4. Varying the FUV Radiation Field 



The FUV radiation field is a significant source of heating of neutral interstellar gas. As 
explained in § (3J the radiation field is taken to be the standard interstellar radiation field 
( lDraindll978l ) multiplied by a factor of \ where x = 1-7 corresponds to a flux of 1.6 x 10~ 3 
erg s _1 cm -2 . The FUV radiation field strength is linked to the number of m assive stars 
in the Galaxy and given the presence of many massive stars at high redshift (lAbel et al. 



20021 ; iBromm et al.lll999l ) is likely to be large in these systems. In this study, we vary x 
between 1.7 and 1.7xl0 4 and study the effect on our star formation timescales. The results 
are summarised in Figure HI 

At A v ~ 10, the FUV radiation is sufficiently attenuated to not affect the chemistry 
within the cloud and all our timescales therefore remain fairly constant with FUV flux 
strength. At lower visual extinctions, we find that the cooling time is almost always smaller 
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than the free-fall time and collapse and star formation are therefore encouraged. The freeze- 
out time is also almost always smaller than the desorption timescale except at very high 
FUV radiation field strengths. This suggests that the formation of icy mantles on dust 
grains is possible even at FUV radiation field strengths of % = 1.7 x 10 3 at A v ~ 3. The 
desorption time and freeze-out time are fairly comparable at all radiation field strengths 
so any icy mantles deposited on the dust grain are quick to evaporate. At very high FUV 
flux strengths, the desorption time is shorter than the freeze-out time and ice formation is 
inhibited. Both these processes mean that the ISM is always rich in molecular coolants at 
high FUV radiation field strengths and collapse is encouraged. 

The ambipolar diffusion timescale is very large for high radiation field strengths and at 
low A v . This once again suggests that magnetically super-critical cores will be stable against 
collapse. However, at high A v where the radiation field cannot penetrate the cloud, the 
ambipolar diffusion timescale is comparable to the dynamical timescale and the magnetical 
criticality is less significant. 

As we are linking both the FUV radiation field strengths and the cosmic ray ionisation 
rate with massive stars in this rather simplistic analysis, it is not surprising that both 
parameters have similar effects on the star formation timescales. We conclude that the 
presence of massive stars resulting in an elevated FUV flux strength, serves to encourage 
quiescent low-mass star formation in every possible way at low A v provided the residual 
magnetic pressure implied by the long ambipolar diffusion timescale, does not halt collapse. 
We compute the ratio of the magnetic to thermal pressure for \ = 1-7 x 10 4 to be a ~ 0.5 
and conclude that the magnetic pressure may in some cases be sufficient to halt collapse. 



4.5. Modelling a High Redshift Object 



Bavet et al.l (120081 ) have identified the Cloverleaf (z ~ 2.6) and APM08279 (z ~ 3.9) as 



possible models for high redshift sources. Although the physical properties of these sources 
are not well constrained, molecular detections suggest th at they are actively forming stars 
(ILutz et al.l 120071 ; iGao &: Solomon! 12004 ; IWagg et al.ll2005l ) . We expect such objects to have 
a range of different metall icities with galaxie s with significant observed dust obscuration 
such as HCM6A (z = 6.56 JChary et al.l (120051 )) likely to be more chemically enriched. The 
presence of many massive stars in these systems will also result in high FUV radiation field 
strengths and cosmic ray ionisation rates. Collapse is likely to occur at moderately high 
densities due to the lack of available molecular coolants compared to a Milky Way type 
object. 
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In order to model such a system, we calculate the star formation timescales for three 
different models - (I) n = 10 4 cm" 3 , f = O.O5£ , C = 10~ 13 s' 1 and x = 1-7 x 10 4 , (II) 
n = 10 5 cm" 3 , £ = O.O5£ , ( = 10" 14 s" 1 and x = 1.7 x 10 3 and (III) n = 10 5 cm" 3 , 
£ — £©; C — 10~ 14 s _1 and ^ = 1.7 x 10 3 - at depths corresponding to A v ~ 3 and A„ ~ 10. 
As previously stated, the presence of many massive stars in high redshift systems justifies 
our choice of elevated FUV radiation field strengths when modelling s uch systems. In order 



to mimic star burst galaxies with hi gh cosmic ray fluxes such as M82 (ISuchkov et al.lll993l ) 



and N GC253 ([Bradford et al.l 120031 ) at high redshift, we adopt high cosmic ray ionisation 



Meijerink et al. estimate a cosmic ray ionization rate of 5 x 10 15 s 1 for a 



rates. 

star formation rate of ~100M Q yr _1 and we expect high redshift galaxies with higher star 
formation rates to have elevated cosmic ray fluxes. Finally, there is evidence that high 
redshift galaxies exist at a range of metallicities with galaxies with significant observed dust 
obscuration likely to be more chemically enriched and we therefore consider both solar and 
sub-solar metallicities. Our results are summarised in Table [TJ 

We can see that for the first model, the free-fall time is comparable to the cooling time 
and cores can collapse to form stars even at depths as low as A v ~ 3. The formation of icy 
mantles is likely to be inhibited in such systems as the freeze-out time is very large compared 
to the desorption time. 

The second model has slightly lower cosmic ray fluxes and FUV radiation field strengths 
but is considered to be at a higher density. In this case, the free-fall timescale is shorter than 
the cooling time at low A v and cores can collapse to form stars but at A v ~ 10, the cooling 
time is larger and the propensity for low mass star formation is reduced. Once again, the 
formation of icy mantles is inhibited in such a system. 

The third model also has elevated FUV radiation field strengths and cosmic ray fluxes 
but is now considered at solar metallicity. The cooling time is now significantly lower than 
the free-fall time at all depths and low-mass star formation is unimpeded by thermal pressure 
while ice formation still remains inhibited. 

In all three cases, the ambipolar diffusion timescale is very large so the residual field 
may be able to halt collapse. Using the temperatures and densities of these model runs, we 
calculate the thermal and magnetic pressure for each of the models in the same way as in 
§ 14.31 and their ratio, a. In both the sub-solar metallicity models, we compute a ~ 0.05 
so the thermal pressure dominates over the magnetic pressure. The cooling times are small 
compared to the dynamical free-fall time, particularly at low A v and low-mass star formation 
is encouraged. The third model at solar metallicity has a ~ 0.6 suggesting that in this case 
the magnetic pressure may be sufficient to halt collapse. We conclude that recent evidence 
for a top-heavy IMF is only supported by this study if high-redshift galaxies are considered 
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to be at solar metallicities in which case the magnetic pressure may halt the collapse and 
formation of low- mass stars. 



5. Conclusions 

In this work we have investigated the variation of certain characteristic timescales as- 
sociated with low mass star formation with physical properties of the collapsing cloud. By 
looking at the balance of these timescales we are able to determine the conditions necessary 
for the formation of low mass stars and comment on the form of the stellar initial mass 
function in extragalactic environments. Despite the crudity of our calculations, we observe 
the following trends: 

(i) Increasing the density results in the cooling time becoming longer than the free-fall 
time, so that cores tend to heat up and star formation is inhibited. Assuming that the 
microscopic properties of dust are similar in extragalactic environments to that in our own 
galaxy, we find that the formation of molecular ices is likely to occur at all densities between 
100 cm -3 and 10 6 cm -3 . The ambipolar diffusion timescale is comparable to the free-fall time 
except at very low densities where it is much longer, suggesting that low mass star-formation 
is allowed in magnetically sub-critical clouds with densities above 10 3 cm" 3 . 

(ii) Low metallicity systems are characterised by long cooling times compared to the 
free-fall timescale. Collapsing cores tend to heat up thus inhibiting star formation and 
collapse. Due to the low dust grain number density, the desorption time is much shorter 
than the freeze-out time in low metallicity systems and there is less ice formation than in 
the Milky Way. 

(iii) High cosmic ray ionisation rates and high FUV radiation field strengths are associ- 
ated with the presence of massive stars. Such systems show rapid cooling times compared to 
their free-fall timescales and star formation is unimpeded by thermal pressure. The propen- 
sity for ice formation is reduced in systems with cosmic ray ionisation rates that are much 
larger than those in the Milky Way. Freeze-out occurs over long timescales and if any ice 
was formed, it would evaporate quickly; molecular coolants would be returned to the ISM 
further increasing the star formation efficiency. The freeze-out and desorption timescales are 
comparable for systems at almost all FUV radiation field strengths so that icy mantles are 
able to form except at very high values of \. The ambipolar diffusion timescale is very long 
for high values of the cosmic ray ionisation rate and radiation field strength. This implies 
that star-formation may not be viable in clouds that are magnetically sub-critical under 
these conditions. 
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(iv) High redshift systems are characterised by high cosmic ray fluxes and FUV radiation 
field strengths and can be considered to exist at a range of metallicities. All such systems 
will show enhanced low-mass star formation with the propensity for star formation increas- 
ing with the metallicity of the galaxy provided thermal pressure dominates the magnetic 
pressure. This is certainly shown to be the case in active high redshift galaxies at sub-solar 
metallicities. However, in active galaxies at solar metallicities, the magnetic field can impede 
the formation of low-mass stars resulting in a high-mass biased stellar IMF. 

While this study is relatively crude in nature, it is instructive in terms of making some 
predictions about essentially unobservable low-mass star formation in extragalactic envi- 
ronments that are different from our own Galaxy. We have noted the kinds of systems in 
which star formation is likely to be enhanced or inhibited by considering changes in various 
timescales relevant for star formation. We found among other things that low-mass star 
formation is likely to be encouraged in host galaxies of massive stars and AGN and inhibited 
in low metallicity systems with low cosmic ray fluxes and FUV radiation field strengths. 
Our study finds tenuous evidence for a top-heavy stellar IMF at high redshifts but only if 
galaxies at high redshift are considered to be at near solar metallicities. In the future, we 
hope to extend this study by varying other free parameters such as the individual elemental 
abundances. We also hope to be able to relate this study to empirical estimations of star 
formation rates in high redshift systems. 

We acknowledge the referee for insightful comments and helping to place the paper 
within a broader context. MB is supported by an STFC studentship. 
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Table 1: Star Formation Timescales at High-Redshift 



Timescale 


Model I 


Model II 


Model III 


Years 


A v ~ 3 A v ~ 10 


A v ~3 A v ~ 10 


A v ~ 3 A v ~ 10 


Free-Fail 
Cooling 
Freeze-out 
Desorption 
Ambipolar Diffusion 
Chemistry 


7.50 x 10 5 7.50 x 10 5 
6.76 x 10 5 7.98 x 10 5 
1.68 x 10 6 1.35 x 10 6 
3.11 x 10 2 1.94 x 10 2 
7.14 x 10 10 8.58 x 10 10 
2.50 x 10 1 2.50 x 10 1 


2.37 x 10 5 2.37 x 10 5 
4.31 x 10 5 6.75 x 10 6 
5.22 x 10 5 1.38 x 10 5 
1.11 x 10 2 2.35 x 10 1 
1.76 x 10 8 2.45 x 10 7 
2.50 x 10 2 2.50 x 10 2 


2.37 x 10 5 2.37 x 10 5 
3.29 x 10 3 1.07 x 10 4 
4.01 x 10 4 3.93 x 10 4 
1.44 x 10 2 1.81 x 10 2 
1.43 x 10 9 1.94 x 10 7 
5.00 x 10 3 5.00 x 10 3 
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Fig. 1. — Variation of timescales relevant to star formation with number density of hydrogen 
at A v ~3 (solid lines) and A v ~10 (dotted lines). The metallicity, cosmic ray ionisation rate 
and incident FUV radiation intensity are fixed at £ = £ Q) ( = 10~ 17 s _1 and \ = 1.7. 
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Fig. 2. — Variation of timescales relevant to star formation with metallicity at A v ~3 (solid 
lines) and A v ~10 (dotted lines). The number density of hydrogen, cosmic ray ionisation rate 
and incident FUV radiation intensity are fixed at n — 10 5 cm -3 , ( = 10 _17 s _1 and \ = 1.7. 
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Fig. 3. — Variation of timescales relevant to star formation with the cosmic ray ionisation 
rate at A v ~3 (solid lines) and A v ~10 (dotted lines). The number density of hydrogen, 
metallicity and incident FUV radiation intensity are fixed at n = 10 5 cm -3 , £ = £o and 
X = l-7. 
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0.5 1 1.5 2 2.5 3 3.5 4 
log 10 (FUV Radiation Field Strength/%) 

Fig. 4. — Variation of timescales relevant to star formation with the FUV radiation field 
strength at A v ~3 (solid lines) and A v ~10 (dotted lines). The number density of hydrogen, 
metallicity and cosmic ray ionisation rate are fixed at n = 10 5 cm~ 3 , £ = £ and ( = 
1 x 10- 17 5 _1 . 
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